ITEP-PH-11/99 



The quark and gluon condensates and low-energy QCD 
theorems in a magnetic field 



N.O. Agasian [] 

and 

0\ ■ I. A. Shushpanov PI 

G\ ■ 

Institute for Theoretical and Experimental Physics, B. Cheremushkinskaya 25, Moscow 

>\ 117218, Russia 

O 

November 1999 

(N 
> 

(N 1 Abstract 



a: 



- 1—1 

X 



The low-energy QCD theorems are generalized in the presence of a constant 
magnetic field H. Two-loop approximation for the vacuum energy density in the 
framework of the chiral perturbation theory was obtained and the quark and gluon 
condensates were found as the functions of H. 



a.: 

• 1 Introduction. 



The low-energy theorems, playing an important role in the understanding the vacuum 



state properties in quantum field theory, were discovered almost at the same as quan- 
tum field methods appeared in particle physics (see, for example, Low theorems |T|). In 
QCD, these theorems were discovered in the beginning of eighties 0. These theorems, 
being derived from the very general symmetrical considerations and not depending on 
the details of confinement mechanism, sometimes gives information which is not easy to 
obtain in another way. Also, they can be used as "physically sensible" restrictions in the 
constructing of effective theories. Recently, they were generalized to finite temperature 
and chemical potential case 0]. 

The behavior of the quark condensate in the presence of a magnetic field was studied 
in Nambu-Jona-Lasinio model earlier For QCD, the analogous investigation in one- 
loop approximations was done in 0. It was found that the quark condensate grows with 
the increase of the magnetic field H in both cases. It implies that naive analogy with 

1 e-mail: agasyan@heron.itep.ru 
2 e-mail: shushpan@heron.itep.ru 



1 



superconductivity, where the order parameter vanishes at some critical field, is not valid 
here. 

In this paper, we generalize the low-energy theorems to the case of the presence of 
the constant magnetic field. When the field is weak (compared to the characteristic 
hadronic scale), two- loop formulae for the vacuum energy density, the quark and gluon 
condensates as the functions of H were obtained basing on the chiral perturbation theory 
(ChPT). Gluons do not carry electric charge, nevertheless, virtual quarks produced by 
them interact with electromagnetic field and lead to the changes in the gluon condensate. 



2 Low-energy theorems in the magnetic field. 

The Euclidean version of QCD partition function in the magnetic field can be written 
as follows 

Z = exp J d 4 xF 2 ^ J [DB;] [Dq f ] [Vq f ] exp{- j d 4 xL QCD }, (1) 
where QCD Lagrangian in the background field is 

\ a 1 

Qf, (2) 



Here Qf is the charge matrix of quarks in the flavor space and we have suppressed gauge 
fixed and Faddeev-Popov terms since they are unessential here. The vacuum energy 
density is defined in the usual way, Vie v (H,mf) = —InZ. In the constant magnetic field, 
the chiral condensate in the chiral limit (m/ — > 0) is given by 



{mi )(H)= eUH ' mi) 



dm* 



(3) 



mt=0 



From (HD follows that the gluon condensate (G 2 ) = {(G^ u ) 2 ) can be written in the following 
form 

( G 2^ = a de v (H,m f ) 



d(l/g 2 ) 



(4) 

m/=0 



The effect of dimensional transmutation produces a new dimensionful nonperturbative 
parameter 

da* 



A = Mexp J r 

[Ja s (M 



(5) 

S (M) P(a 8 ) J ' 

where M - is ultraviolet cut-off, a s = g 2 /^, and f3(a s ) = da s (M)/d\nM is Gell-Mann - 
Low function. 

The system, corresponding to partition function ([!]), contains two dimensionful param- 
eters M and H in the chiral limit (rrif = 0). As the vacuum energy density is an observable 
quantity, it is also renorm-invariant one and its anomalous dimension is zero. Therefore, 
e v has only normal (canonical) dimension which equals four. Using the renorm-invariance 
of A, we can write the most general expression for e v 

e v = A 4 f(H/A 2 ), (6) 
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where f is some function. Then, it is easy to get from Eqs. (|5]) and that 

d{l/9l) ~ %) { dH) e ^ {7) 
and from Eq. (|), we can find the relation between e v and the gluon condensate 

If we set H = we would reproduce the well-known expression for nonperturbative density 
of the vacuum energy in the chiral limit. In one-loop approximation {(3 = — 6o«s/2vr, 
6 = (IIN C - 2N f )/3), it reads 

b o mi 



128tt ; 



m- (9) 



Let us now derive the low-energy theorems in the presence of magnetic field. We can 
iterate the procedure described above by taking n further derivatives of Eq. (fD 

(_ ir ( 4 - 2ff _j e„=( 2 //--4) („> = 

d 4 x n ■ ■ ■ d 4 X\ (cr(x n ) ■ ■ -a(xi)) c . (10) 

Here 

a(x) = 9 m {x) = Y^- 2 (G%(x)) 2 , (11) 



and the subscript 'c' means that only connected diagrams are to be included. Proceeding 
in the same way, it is possible to obtain the similar theorems for renorm-invariant operator 
0(x) which is built from quark and/or gluon fields 



^2H^-dj (0)= Jd 4 x n ---d l 



x 1 (a(x n )---a(x 1 )O{0)) c , (12) 



where d - is canonical dimension of operator O. If operator O has also anomalous dimen- 
sion, the appropriate 7-function should be included. 



3 Vacuum energy in the magnetic field. 

The obtained results permit us to calculate the condensates in the chiral limit as functions 
of H, once the vacuum energy density is known. To calculate e v , we have to consider the 
loops in the magnetic field. When the field is weak, eH << [i\ adr ~ (AtcF^) 2 , characteristic 
momenta in loops are small and theory is adequately described by the low-energy effective 
Lagrangian |J which admits an expansion in powers of the external momenta (derivatives) 
and the quark masses [] 

L e// = L< 2 ) + I (4) +L (6) + --- (13) 
3 Note that the chiral limit means eH >> M"i. 
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The leading term in (|I5D is similar to the Lagrangian of nonlinear cr-model in external 



field V„ 



L (2) 



4 



-Tr(V M [/ + V M C/) + EReT^jMCT 



(14) 



Here C/ stands for an unitary SU (2) matrix, F n = 93 Mev is the pion decay constant and 
parameter E has the meaning of the quark condensate £ = | (mm) | = | (dd) \ . The external 
Abelian magnetic field H, aligned along z-axis corresponds to V^(x) = e(r 3 /2)y4 M (x), 
where vector potential may be chosen as Ai(x) = —Hx2- 

The difference of light quark mass m u — enters in the effective Lagrangian fll~3|) only 
quadratically. To calculate the fermion condensate, we need to take only first derivative 
over the quark mass before going to the chiral limit. In the case of the gluon condensate, 
we can put m u = irtd = from the very beginning. It means that we can take mass matrix 
to be diagonal M. = ml. 



In the leading approximation the Lagrangian ( |T4D coincides with the Lagrangian of 
scalar electrodynamics. In one-loop approximation, the result for the vacuum energy 
density in this theory was found long time ago by Schwinger fTj 



(H) 



1 



16vr 2 



ds 



-Mis 



eHs 



sinh eHs 



- 1 



(15) 



Using Eqs. (§), (|) and Gell-Mann - Oakes - Renner relation (GMOR) 2mS = F 2 M 2 , we 
can get one-loop formulae for the condensates in the constant magnetic field 



E(H) = E 



1 



e#ln2 
(4^F 



(16) 



(G 2 )(H) = (G 2 ) + 



at 



(eHf. 



(17) 



3nf3(a s 

The expression ([H]) for the quark condensate was obtained in || earlier. In the chiral 
perturbation theory in the magnetic field, the parameter of the expansion is eHj (47rF v ) 2 . 
To find e v to the next order or the expansion, we have to take into account two-loop 
diagrams with the vertices from L^ 2 \ one-loop diagrams with one vertex from and 
a tree contribution arising from L^ 6 \ The corresponded Feynman graphs for e^, 2 -* are 
depicted in Fig. 1. 

To perform the calculations of these diagrams, we should expand up to four- 
pion vertices. The matrix U can be parameterized in many ways. We choose Weinberg 
parameterization 



U = a + 



7T 

1 7T 



1. 



;is) 



Then the expansion of is 



M 2 (tt ) 2 



+ {d^ + + ieA M 7r + )((9^7r - ieA^~ 



M 2 



7T 7T 



4 



/ 

I 

\ 

\ 




(a) 



00 





(c) 



(d) 



Figure 1: The loop diagrams contribute to vacuum energy density to the second order 
of ChPT. Solid lines correspond to the charged pions and dashed lines correspond to the 
neutral pions. 



1 M 2 
+ ^[n°dy + ^(tt+tt-)] 2 - ^[27r+7r- 



7T 



0\2i2 



2F 2i n 8F 2< 

where we have introduced the fields of the charged and the neutral pions 

1 



7T° = 7T 3 , 



7T~ 



V2 



(tt 1 ±in 2 ). 



(19) 



(20) 



We need also the expression of a charged pion propagator in a magnetic field. It can 
be inferred from the results of Ref. || where an explicit expression for the fermion 
propagator in a magnetic field at nonzero chemical potential \i has been found. After some 
simplifications, we arrive at the following formula for the Euclidean scalar propagator 



D H (x,y) = *(z,y) J -0^e lk ^D H (k) 



(21) 



where $(x, y) = exp{ie Jy A lx {z)dz IJ } is the phase factor and the integration is done along 
the straight line connecting x and y, and 



D H {k) 



ds 



o cosh(eHs) 



i / l9 , 9 tanheifs , 
exp {—s\k?\+k A + Ml 



eHs 



(22) 



with k 2 = k\ + kj, k\ = kj + k\. 

The contribution of the graph in Fig. la to e v is 



4 2) [Fig-la] 



M. 



^ 2 (o), 



(23) 
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where D(0) is a free propagator of the scalar particle at coinciding points which can be 
written in dimensional regularization scheme as 

f d d k 1 „ / 1 M 2 \ , . 

where \x is the mass parameter of the regularization and a singular term in (£23|) is ex- 
plicitely separated out as 



A 



a d - A r 1 1 , 

■(ln47r -7e + 1) 



16vr 2 



d-4 2 



(25) 



The correction (^3|) can be absorbed to e v at = and does not cause the shifts of the 
condensates by magnetic field. 

The correction to e v , coming from the next diagram Fig. lb, is 

M 2 

^Wb] = ^£(0)^(0), (26) 

where D H (0) = D (x, x). According to GMOR, the result in (|26|) is proportional to the 
square of the quark mass and does not change the condensates in the chiral limit. 
Performing the calculation of the diagram Fig. lc, we get 

e^Fig.lc] = j- 2 D H (0) J ^(k 2 + M 2 )D H (k). (27) 

This expression contains a quartic divergence. However, by virtue of the dimensional 
regularization identity J d d k = 0, this divergence can be ignored 0. Subtracting one from 
the integrand in ( |27| ) and going to a limit d — > 4, we find that [Fig.lc] = 0. So, we 
arrive at the conclusion that two-loop diagrams do not give the correction to e v , linear in 
the quark mass and, therefore, do not shift the condensates. 

In this order of expansion in eH/(47rF n ) 2 , there are some diagrams with vertices from 
besides the diagrams discussed above. As the momentum of a constant field is zero, 
only the two terms from the general form of are left 

= -^(ei^vr+TT- - ^eF^n-d^ + ieA,d v (ir+7r-)}, (28) 

7T 7T 

where the phenomenological constants (infinite) I5 Iq were defined in 0. The correspond- 
ing diagrams are drawn in Fig. Id. The calculation is straightforward and gives 

4 2) [Fig.ld] = ^£(2/ 5 -/ 6 )^(0). (29) 

Although the constants Zg l@ are infinite, their combination, arising in (|29|) is finite B 



2l 5 -k = -^(k-l 5 ) J (30) 



4 For other regularizations, quartic divergence is removed only if a nontrivial measure for the group 
integration is taken into account. The examples of explicit calculations can be found in H. 
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and l 6 — h ~ 2.7. 

Total expression given by chiral symmetry, is rather complicated ||10|| . However, 
for our purposes, the only one term is important which can be taken in the following form, 
convenient for us 



L (6) 



80d 



{eF llu ) 2 YReTi{MU + } 



(31) 



Here d = d r {jj) + const ■ A and A involves a pole ^ d ~^ /{d — 4). The divergent part in (|3l"D 
is canceled by the pole terms, stemming from one-loop diagrams with vertices from . 
Numerical value of d r can be inferred from the results obtained in [1 i| where the process 
77 — > 7T 7r° has been studied. In the notations of d r can be rewritten as 
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d r 



320 



(ai + 2a 2 + 46), 



(32) 



where a\, d2 and h have the meaning of the coefficients at different tensor and mass 
structures in the amplitude of the process 77 — > 7r°7r°. The numerical values of ax, a<i and 
b have been obtained by vector, tensor and scalar resonances saturation of the amplitude. 



11]. It turns out that only the exchange by scalar mesons gives the contribution to d r 
and 

cf Gu ~ 0.75 Gev) « ±4 ■ 10~ 6 (33) 

The scalar meson-photons and meson-pions coupling constants enter in the experimentally 
observable decay widths quadratically, since in the effective chiral Lagrangian they emerge 
linearly. Hence, it easy to see that the sign of d r remains undetermined. 

Now we are in the position to get a finite result for e v . To find the quark condensate, 
one should keep only first order terms in M 2 / eH expansion 



^(0) 
Ml 



[D H (0)-D(0)] + D(0) 



eg In 2 



+ 



16vr 2 



1 e H „ 



+ 2Mi 



167T 2 

1 M 2 
A + -— ln-f 



32tt 



(34) 



where C - is slowly varying function of eif/M 2 and C(0) ~ —0.2. Adding up the obtained 
results, we arrive at the final answer for e v in two-loop approximation 



(eHf 



48vr 2 (4vrF^) 2 
160d r (fi) 



(h-h)ly-eH\n2 + Ml 
eHfMl. 



1 e H „ 
In — + C 

fi 2 



(35) 



Using (^ 
upon H 



and GMOR relation, it is easy to find how the quark condensate depends 



£(#) = S 1 + 



eH 



■In 2 - - 



1 {eHf 



(4tt^) 2 " 3(4ttF t ) 4 



(r 6 




(36) 
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Figure 2: The dependence of the gluon condensate (G 2 (x)) / (G 2 ) (solid line) and the 
quark condensate (dotted line corresponds to d r > and dashed line corresponds 

to d r < 0) upon eH/(AnF w ) 2 . 



Further, we have for the gluon condensate 



(G 2 )(H) = (G 2 ) + 



3ir(3(a s 



{eHf 



(37) 



It follows from asymptotic freedom that (3{a s ) < and, consequently, the gluon conden- 
sate drops when the magnetic field arises. In one loop approximation for QCD (3- function 
(/3(a s ) = — boa 2 /2ii ), this decrease is equal to 



A(G 2 ) 



2vr 
3b~o 



eHf 



1 + 2 



eH 



(4vriV 



-(r 6 -r 5 ) in 2 



(38) 



Let us introduce the dimensionless variable x = eH / [AtxF^) 2 . Substituting the found 
numerical values into Eq. (|3B|), we can write the quark condensate as a function of x 

£(z)/£ = l + x\n2-ax 2 \nx-bx 2 . (39) 

Here a ~ 0.9 and b ~ 0.65 ± 1.77, where, for the coefficient b, we have used two values of 
d r from (p3|). The behavior of condensates (G 2 (x))/(G 2 ) and S(x)/E (for two different 
values of b) is shown in Fig. 2. If d r > then the quark condensate starts to decrease at 
the magnitude of the magnetic field x > 0.23. 



4 Conclusions 

We have generalized the low-energy QCD theorems in the presence of the constant mag- 
netic field H. Basing on ChPT, we have calculated the vacuum energy density in two-loop 
approximation and found the dependence of the quark and gluon condensates upon the 
magnitude of H. It is shown that the gluon condensate diminishes with the increasing 
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of H while the behavior of the quark condensate crucially depends on the sign of d r . It 
decreases if we choose positive sign and increases in the opposite case. Note that pos- 
sible decreasing of the condensate S happens in the region of the applicability of ChPT 
eH/^TrFn) 2 < 1. As it was already mentioned, it is not possible to determine the sign 
of d r from experimental data. We will not bring here various speculations concerning the 
behavior of £ in the magnetic field. We would like only to say that the question about the 
sign of E changing in two-loop level would be resolved if the phenomenological constants 
of ChPT are derived from the first principles of QCD. 
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